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SUMMARY 

The>  development  of  our  knowledge  of  the  bmckling  of  thin-walled 
circular  cylindrical  shells  subjected  to  axial  compression  is  outlined5 
from  the  beginning  of  the  century  until  the  present,  with  particular 
emphasis  on  advances  made  in  the  last  twenty -five  years.  It  is  shown 
that  practical  shells  generally  buckle  under  stresses  much  smaller  than 
the  classical  critical  value  derived  by  Lorenz,  Timoshenko,  Southwell 
and  Flugge.  A  first  explanation  of  the  reasons  for  the  discrepancy  was 
given  by  Donnell  and  the  problem  was  explored  in  detail  by  von  Karman, 
Tsien  and  their  collaborators.  More  recently,  Yoshimura  discovered  the 
existence  of  an  inextensional  displacement  pattern  which  the  wall  of 
the  shell  can  suddenly  assume,  and  Koiter  found  an  explanation  of  the 
sensitivity  of  the  buckling  stress  to<  small  initial  deviations  from  the 
exact  circular  cylindrical  shape. 

In  the  last  few  years  further  interesting  discoveries  were  made 
in  Japan  and  in  California  regarding  the  effects  of  details  of  the 
boundary  conditions,  and  many  additional  numerical  results  were  obtained 
with  the  aid  of  high-speed  electronic  digital  computers.  Improvements 
in  experimental  techniques  have  also  contributed  significantly  to  a 
clarification  of  the  problem  and  to  an  establishment  of  the  unavoidable 
deviations  from  the  exact  shape  as  the  major  causes  of  the  .large  differ¬ 
ences  between  theory  and  experiment. 
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INTRODUCTION* 

Theodore  von-  Karman,  in  whose  memory  this  lecture  is  "being  pre¬ 
sented,  is  best  known  for  his  work  in  aerodynamics.  Nevertheless  it 
is  appropriate  to  speak  about  structures  in  a  von  ;Karraan  Memorial 
Lecture  because  few,  if  any,  research  men  have  contributed  as  much  as 
did  von  Karman  to  the  solution  of  problems  of  structural  stability. 

His  doctoral  thesis  (von  Karman  1908,  1910a)  is  the  foundation  of  the 
theory  of  buckling  at  stresses  exceeding  the  elastic  limit  of  the 
material.  It  shows  that  the  buckling  load  of  a  short  column  can  be 
calculated  from  a  modified  Euler  formula  in  which  Young’s  modulus  of 
elasticity  is  replaced  by  a  combination  of  the  tangent  modulus  and 
Young’s  modulus,  generally  known  as  the  von  Karman  modulus.  In  spite- 
of  some  recent  attempts  to  replace  the  von  Karman  buckling  load  with 
the  tangent -modulus  load,  it  is  now  perfectly  clear  (Hoff  1965d)  t-hat- 
the  von  Karman  load  is,  and  will  always  be,  the  correct  and  the  only 
solution  of  the  buckling  problem  of  short  columns  in  the  classical, 
or  Eulerian,  sense. 

An  equally  important,  and  even  more  original,  contribution  of  the 
man  whom  we  honor  today  is  the  explanation  of  the  snap-through  type  of 
buckling  of  thin  shells  subjected  to  compression.  In  a  series  of  papers 
(von  Karman  et  al.  1939>  19^0,  19^1)  he  an  his  collaborators  proved  the 
existence  of  equilibrium  states  corresponding  to  large  deformations  of 
the  compressed  shell.  These  states  are  responsible  for  the  low  experi¬ 
mental  vaLues  obtained  and  for  the  perplexing  behavior  of  the  shell 
when  it  buckles.  These  topics  will  be  discussed  at  some  length  in  the 
body  of  the  paper. 
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But  it  would  be  foolish  to  try  to  assess  the  importance  of  von 
Harman' s  work  on  the  basis  of  his  published  contributions  to  aerodynamics 
and  structural  stability.  The  111  original  papers  reprinted  on  the  1838 
pages  of  his  -Collected'  Works  (von  Karman  1956)  deal  also  wish  stress 
analysis,  plasticity,  hydrodynamics,  flight  mechanics,  jet  and.  rocket 
propulsion,  combustion,  physics,  mathematics,  and  airplane,  missile  and 
machine  design.  They  were  written  in  four  languages,  English,  German, 
French  and  Hungarian.  The  mother  tongues  of  Dr.  von  Karman1 s  friends, 
collaborators  and  students  are  even  more  numerous,  and  the  influence  of 
von  Karman  on  mechanics,  aeronautics  and  astronautics  will  always  be 
felt  in  all  the  countries  of  the  world. 
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THE  .DEVELOPMENT  OF  THE  CLASSICAL  BUCKLING  FORMULA 


kling 

The  classical  "buckling  formula  is  usually  written  in  the  form 
--l/2 

°cl  =  &(1'V2  J  E(h/a)  «  0.6s(h/a)  (l) 

where  c  .  is  the  critical  value  of  the  uniform  axial  stress.  E  is 
cl  f 

Young's  modulus  of  elasticity ,  v  Poisson's  ratio,  h  the  thickness 
of  the  wall  of  the  circular  cylindrical  shell,  and  a  the  radius  of 
the  middle  surface  of  the  shell.  Essentially  this  expression  was 
derived  by  Rudolf  Lorenz,  a  civil  engineer  in  Dortmund,  Germany,  in 
1908.  Starting  out  from  the  expressions  in  volume  5  of  Foppl's  Tech- 
nische  Mechanik  (Foppl  1898),  Lorenz  developed  the  equations  governing 
the  axisymmetric  deformations  of  initially  .slightly  inaccurate  thin- 
walled  circular  cylindrical  shells  subjected  to  uniform  axial  compression. 
For  simply  supported  circular  edges  he  represented  both  the  initial 
deviations  from  the  exact  shape  and  the  additional  displacements  due  to 
the  load  by  Fourier  sine  series  and  obtained  the  critical  value  of  the 
stress  as  the  one  at  which  the  denominator  of  one  of  the  particular 
solutions  vanished,  thus  Implying  an  increase  of  the  deformations  beyond 
all  bounds.  He  also  found  a  good  approximation  to  the  correct  expression 


for  the  wave  length 


X  =  Tr|l2(l-v2)]  ^  (ah)1/2  a  l*72(ah)1'/2 


Since  this  wave  length  is  small  for  thin-walled  shells,  and  since 
a  small  change  in  the  relatively  large  number  of  waves  in  the  axial 
direction  changes  but  slightly  the  buckling  stress,  Eqs.(l)  and  (2)  can 


be  assumed  to  represent  with  engineering  accuracy  the  buckling  stress 
and  the  wave  length  of  moderately  long  shells.  The  expressions  become 
rigorously  correct  in  the  limit  when  the  length  L*  increases  indefi¬ 
nitely,  but,  of  course,  very  long  shells  can  buckle  as  Euler  columns, 
without  any  bulges  or  waves  developing  in.  the  shell  wall. 

Lorenz1  formulas  can  be  obtained  from  the  two  equations  given  here 
by  setting  Poisson's  ratio  equal  to  zero.  Thus  his  buckling  stress  was 
too  low  by  5  percent,  and  his  wave  length  was  too  short  by  2§  percent. 
The  correct  formulas  were  first  given  in  a  Western  European  language  by 
Timoshenko  in  1910,  although  they  were  probably  published  earlier  in 
Russian.  In  his  Zeitschrift  fur  Mathematik  und  Physik  article  the 
buckling  stress  is  derived  once  by  the  energy  method,  and-  a  second  time 
through  solution  of  an  eigenvalue  problem  defined  by  a  fourth-order 
ordinary  differential  equation  and  suitable  boundary  conditions,  and 
identical  results  are  obtained  by  the  two  approaches. 

Chessboard  Type  of  Buckling 

The  general  case  of  buckling,  without  the  restriction  to  axially 
symmetric  deformations,  was  first  tackled  by  Rudolf  Lorenz  in  1911.  He 
started  out  from  Love's  so-called  first  approximation  shell  theory  and 
assumed  that  the  displacements  in  the  axial  direction,  designated  as 
the  u*  displacements  (see  Fig,  l),  were  negligibly  small.  He  arrived 
at  a  sixth-order  partial  differential  equation  which  he  solved  in  the 
presence  of  the  following  boundary  conditions: 

v*  =  w*  =  M  =0  when  x*  =  0,L*  (3) 

X 
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Here  v*  and  w*  are  the  displacements  in  the  circumferential  and 
radial  directions,  is  the  axial  "bending  moment  resultant-,  x*  the 

axial  coordinate,  and  L*  the  length  of  the  shell.  He  also  established 
that  the  inclusion  of  the  u*  displacement  in  the  theory  would  yield  an 
eighth-order  differential  equation  whose  solution  would  not  differ 
appreciably  from  that  of  the  sixth-order  equation. 

The  numerical  results  of  the  calculations  were  presented  in  graph¬ 
ical  form.  They  are  in  good  agreement  with  values  computed  from  Eq. (l) . 

A  rather  complete  treatment  of  the  problem  of  elastic  stability  was 
given  by  Southwell  in  1914.  Starting  from  Love's  theory  of  thin  shells, 
he  derived  and  solved  the  equations  defining  the  neutral  equilibrium  of 
a  thin  circular  cylindrical  shell  subjected  to  simultaneous  axial  com¬ 
pression  and  lateral  pressure.  The  same  problem  was  treated  again  by 
Flugge  in  1932  without  reference  to  Love's  theory.  It  follows  from 
both  treatises  that  Eqs.(l)  and  (2)  define  the  buckling  stress  of  the 
axially  compressed  shell  of  moderate  length  in  good  approximation,  and 
that  of  the  infinitely  long  shell  accurately,  provided  that  buckling 
takes  place  symmetrically  to  the  axis  (see  Fig.  2).  When  buckling  is 
of  the  chessboard  type  (Fig.  3),  Eq.(l)  is  still  correct  with  the 
limitations  given,  but  Eq. (2)  is  replaced  by  a  condition  connecting 
the  axial  wave  length  with  the  circumferential  wave  length.  The  wave 
lengths  are  thus  indeterminate;  all  that  can  be  said  is  that  knowledge 
of  one  wave  length  allows  the  calculation  of  the  other. 

These  conclusions,  and  a  comparison  of  the  theory  with  test  results, 
can  be  found  in  Timoshenko's  Theory  of  Elastic  Stability  (Timoshenko  1936). 
The  validity  of  Eq.(l)  in  the  case  of  the  chessboard  type  of  buckling  was 
probably  first  established  by  Timoshenko  in  19-l4. 
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CRITICISM  OF  THE  CLASSICAL  RESULTS 


Buckling  Patterns 

Figures  2  and-  3  represent  the  two-  classical  "buckling  patterns 
derived  from  the  linear  theory.  The  former  can  "be  observed  easily  if 
tests  are  made  with  tubes  of  comparatively  small  a/h  ratio,  ,say  30, 
and  if  the  material  of  the  tube  is  capable  of  large  plastic  deformations, 
which  is  the  case,  for  instance,,  with  mild  steel.  The  deformations  shown, 
in-  Fig.  2  are  permanent;  the  pattern  cannot  be  observed  when  the  deform 
mations  are  entirely  elastic.  But  the  theory  of  buckling  presented 
assumes  perfect  elasticity  of  the  material  which  constitutes  an  incon¬ 
sistency  not  easy  to  resolve. 

The  situation  is  even  worse  with  the  chessboard  pattern  shown  in 
Fig.  3,  It  is  the  reproduction  of  a  drawing  and  not  a  photograph  of  an 
actual  test  specimen,  as  is  Fig-.  2.  The  reason  for  not  showing  a  photo¬ 
graph  is  that  the  pattern,  although  clearly  defined  by  the  solution  of 
the  equations  governing  buckling,  has  never  been  observed  in  actual 
experiment.  The  shape  of  a  thin -walled  perfectly  elastic  specimen  after 
buckling  can  be  seen  in  the  photograph  of  Fig,  k.  As  this  differs  con¬ 
siderably  from  the  former  two  patterns,  the  conclusion  must  be  reached 
that  the  agreement  between  buckled  shapes  predicted  by  theory  and  those 
observed  m  experiment  is  very  poor. 

The  Work  of  Flugge 

In  his  classical  Habilitationsschrift,  that  is  formal  lecture  pre¬ 
sented  when  he  was  appointed  Privatdozent*  in  Gottingen,  Flugge  (1932) 

*  Adjunct  Professor 
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compared  the  results  of  his  analysis  with  data  on  experiments  described 
in  the  technical  literature.  As  he  was  unable  to  find  test  results  oh 
axially  compressed  .shells,  he  manufactured  and  tested  a  number  of  rubber 
and  celluloid  cylinders  in  the  Institute  for  Applied  Mechanics  in  the 
University  of  Gottingen.  The  length -to -radius  ratio  of  the  specimens 
varied  between  l.'jS  and  5,  and  their  radius -to-wall  thickness  ratio  from 
90  to  138.  All  the  specimens  buckled  elastically,  and  the  ratio  of  the 
experimental  buckling  stress  to  the  buckling  stress  according  to  Eq. (l) 
ranged  from  0.52  to  0,65  for  the  celluloid  cylinders. 

To  explore  the  possible  causes  of  this  discrepancy,  Flugge  studied 
carefully  the  boundary  conditions  and  the  effects  of  small  initial 
deviations  from  the  exact  shape.  In  the  classical  approach  the  stability 
of  the  equilibrium  is  examined  when  a  perfectly  cylindrical  shell  is  under 
the  effect  of  a  uniformly  distributed  axial  compressive  stress  a  .  But 
if  the  simple  support  boundary  conditions  exist  at  the  moment  of  buckling, 
they  must  have  been  in  existence  during  the  entire  loading  process  during 
which  the  average  axial  stress  was  brought  from  zero  to  its  value  o  . 

Now  it  is  well  known  that  the  shortening  of  the  shell  L*o/E  is  accom¬ 
panied  by  an  increase-  in  the  radius  amounting  to  vao/E  ,  and  evidently 
this  increase  is  prevented  by  the  simple  support  provided  in  the  two  end 
sections  of  the  shell.  Under  an  axial  compressive  stress  0  the  shell 
must  therefore  have  the  shape  shown  in  Fig.  5  before  it  buckles. 

In  his  investigation,  Flugge  solved  the  sixth-order  ordinary  differ¬ 
ential  equation  defining  the  axially  symmetric  deformations  of  a  shell 
whose  end  sections  are  prevented  from  expanding  during  the  loading  process. 
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He  found  that  at  the  beginning  of  the  loading  the  major  portion  of  the 
■shell  is  cylindrical,,  and  the  meridian  is  slightly  curved  only  near  the 
supports.  As  the  intensity  of  loading  approaches  the  critical  value, 
the  curved  region  extends  toward  tne  middle  of  the  cylinder,  several 
waves  appear  along  the  meridian,  and  the  amplitudes  of  the  waves  increase 
Until  they  become  infinitely  large  at  the  classical  critical  value  of 
the  compressive  stress.  Considerably  below  this  value,  however,  the 
stresses  at  the  -crests  of  the  waves  can  be  large  enough  to  cause  perma¬ 
nent  deformations  even  though  the  average  stress  o ,  is  below  the 

elastic  limit  of  the  material. 

Similarly,  Flugge  showed  that  a  small  initial  deviation  from  the 
exact  shape  in  the  form  of  a  sine  function  of  the  axial  coordinate 
multiplied  by  a  sine  function  of  the  circumferential  coordinate  is 
increased  during  loading  and  tends  to  infinitely  large  values  as  the 
classical  critical  stress  is  approached.  From  this  fact  again  Large 
stresses  and  inelastic  deformations  follow  below  the  critical  load. 

Since  Flugge' s  experimental  buckling  loads  differed  relatively 
little  from  the  theoretical  critical  loads,  the  two  studies  just  sketched 
appeared  to  suffice  for  an  explanation  of  the  discrepancy.  But  in  exper¬ 
iments  carried  out  in  connection  with  the  rapid  development  of  thin- 
walled  aluminum  alloy  airplane  structures  in  the  early  nineteen  thirties 
shells  having  radius -to-thickness  ratios  up  to  1500  were  tested  (Lundquist 
1933,  Donnell  1934).  As  these  specimens  often  failed  at  stresses  as  3ow 
as  15  percent  of  the  classical  theoretical  value,  Flugge !s  explanation 
of  the  causes  of  the  disagreement  was  no  longer  sufficient. 
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JEW  SOLUTIONS  OF  THE  CLASSICAL  EQUATIONS 

The  Semi-Infinite  Shell  with,  a  Free  Edge 

In  19 59  the  present  author  read  a  paper  by  Nachbar  in  which  the 
effect  of  pressurization  upon  the  influence  coefficients  of  rotationally 
symmetric  shells  was  evaluated-.  When  applied  to  the  spherical  shell  of 
not  too  small  solid  angle  <p  at  the  edge,  the  result  was  that  the 
influence  coefficients  k. .  ,  that  is  the  generalized  edge  displacements 
caused  by  unit  edge  stress  resultants,  could  be  expressed  as 


k. 

ij 


k 

i  jo  2x  -f- 1 


k, ,  =  k. 


iJ 


ijo  2x  +1 


if  i  =  j 


if  i  4  j 


(ba) 

(4b) 


In  these  equations  k.  is  an  influence  coefficient  in  the  presence  of 
pressurization,  and  k„Q  the  same  influence  coefficient  when  the  shell 
is  not  pressurized o  The  symbol  x  was  defined  in  the  Nomenclature  as 
the  pressurization  parameter  and  was  given  in  terms  of  other  symbols 
defined  elsewhere. 


Since  the  derivations  contained  no  step  that  would  invalidate  the 
results  if  the  internal  pressure  acting  upon  the  shell  were  exchanged 
for  an  external  pressure,  it  was  reasonable  to  assume  that  x  =  -  l/2 
would  correspond  to  the  classical  buckling  stress  given  in  Eq.(l).  This 
equation  is  known  to  be  valid  for  spherical  shells  under  external  pressure, 
and  the  edge  displacements  can  be  expected  to  increase  without  bounds  when 
unit  stress  resultants  are  applied  to  the  edge  in  the  presence  of  the 


s 

.i 
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critical  stress  in  the  shell.  With  x  =  -  l/2  naturally  all  the 
of  Eqs.(ii)  tend  to  infinity. 

When  this  conjecture  was  checked,  it  turned  out  that  x  =  -  1 
rather  than  x  »  -  l/2  corresponded  to  the  classical  critical  value  of 
the  external  pressure.  This  implied,  of  course,  that  the  free  edge  ,of 
a  spherical  cap  would  become  unstable  at  one -half  the  critical  pressure 
of  the  complete  spherical  shell. 

Because  of  the  nature  of  the  governing  equations,  the  same  result 
could  be  anticipated  to  hold  also  in  the  case  of  the  thin-walled  cir¬ 
cular  cylindrical  shell.  A.  simple  derivation  (Hoff  1961)  presented  at 
a  symposium  honoring  Dr.  von  Karman  on  his  80th  anniversary  proved 
indeed  that  the  buckling  stress  of  a  semi-infinite  circular  cylindrical 
shell  whose  near  edge  is  perfectly  free  to  deform  is  one-half  the  classical 
value.  A  unit  length  of  the  free  edge  is  assumed  to  be  subjected  to  a 
st.rest.  resultant  a  h  whose  magnitude  and  direction  remain  unchanged 
during  the  buckling  process  (see  Fig.  6)  in  agreement  with  the  Eulerian 
concept  of  buckling. 

The  eigenvalue  problem  was  stated  in  the  following  form: 

wIV,  (0  /E)4K*w"+  4K4w  =  0  (5) 

where  the  differentiations  indicated  by  superscripts  must  be  carried  out 
with  respect  to  a  non-dimensionalized  axial  coordinate  x  .  This,  and 
the  non-dimensionalized  radial  displacement  w  are  defined  by 

x  =  x*/a  w  =  w*/a  (6) 


-  10  - 


The  following  constants  were,  introduced  in  the  analysis: 


kK *  =  12(1-V2)  (a/h)2  =  h(E/acl)2  D  =  Eh3/l2(l-v2),  (7) 

The  critical  value  o  of  the  uniform  axial  stress  a  is  the  one 

cr  x 

that  satisfies  the  homogeneous  boundary  conditions 


K 

1! 

O 

11 

when 

0 

11 

X 

X 

V  =  -  (D/a2)w'"=  a„  hw' 

'  '  '  cr 

when 

x  =  0 

and  yields  hounded  values  for  all  displacements  and  stresses  when 
x  -*■  «>  in  the  solution  of  Eq.(5).  This  value  was  found  to  he 

o  /o  ,  =  1/2  (9) 

for  the  semi-infinite  shell  that  buckles  axisymmetrically.  It  is  perhaps 

worth  noting  that  the  second  of  Eqs.(8)  expresses  the  fact  that  a  small 

transverse  shear  force  a  h  sin  tan  ^w*  «  .0  hw*  must  act  on  the  edge 

cr  cr 

of  the  shell  if  the  direction  of  the  applied  axial  compressive  load 
remains  unchanged  while  the  generator  at  the  edge  of  the  shell  rotates 
during  buckling. 

The  restriction  of  axisymmetry  in  the  deformations  was  removed  in 
a  follow-up  paper  by  Nachbar  and  Hoff  (Nachbar  et  al.  1961,  19&2)  and 
the  critical  stress  of  the  semi-infinite  shell  was  found  to  be 

0  /o  ,  =  0.38  (10) 

cr  c.1  '  ' 


-  11  - 


Ohira 's  Solution 

In  the  discussion  of  a  lecture  the  author  gave  to  the  Japan  Society 
of  Aeronautics  at  Meiji  University  In  Tokyo  in  April  1963,  Hiroichi  5hira, 
a  young  associate  professor  at  Kyushu  University,  disclosed  that  he  had 
obtained  approximately  l/2  for  p  for  semi-infinite  shells  with  the 
near  edge  simply  supported:  (Ohira  1961).  He  had  discovered  that  changes 
in  the  classical  boundary  conditions  relating  to  displacements  and  stresses 
in  the  tangent  plane  to  the  middle  surface  can  lead  tb  such  a  reduction 
in  the  buckling  stress.  Ohira  used  a  relatively  complicated  differential 
equation  and  obtained  solutions  with  the  aid  of  a  digital  computer. 

The  Donnell  Equations 

It  occurred  to  the  author  of  the  present  paper  that  a  closed-form 
solution  of  the  equations  defining  the  buckling  of  thin-walled  circular 
cylindrical  shells  should  be  possible,  and  that  it  should  yield  the 
same  low  buckling  stress  as  that  obtained  by  Ohira,  if  the  latter's 
solution  was  correct.  To  check  this  conjecture,  he  and  his  graduate 
student  Rehfield  set  up  the  problem  (Hoff  et  al.  1964a,  1965c)  with  the 
aid  of  Donnell's  small  displacement  equations.  These  equations  were 
originally  published  by  Donnell  in  1933,  when  he  was  a  research  associate 
in  the  Guggenheim  Aeronautic  Laboratory  of  the  California  Institute  of 
Technology  whose  director  was  von  Karman.  The  equations  became  popular 
after  Batdorf  had  used  them  extensively  in  the  solution  of  a  number  of 
stability  problems,  and  checked  their  accuracy  against  other  available 
theoretical  solutions  and  test  data  (Batdorf  1947)*  Another  check  of 
their  accuracy  was  presented  by  the  author  when  he  compared  the  roots 


of  the  characteristic  equation  corresponding  to  Donnell’s  equation,  with 
those  obtainable  from  other  shell  equations  available  in  the  literature 
(4off  1955). 


The  Donnell  equations  differ  from  the  many  other  equations  defining 
the  deformations  of  thin  circular  cylindrical  shells,  in  so  far  as  they 
are  much  simpler  and  more  symmetric  in  their  structure  than  most  of  the 
others.  At  the  same  time  they  are  sufficiently  ac-  urate  in  problems  of 
buckling  if  the  number  N  of  waves  around  the  circumference  is  large 
enough.  It  is  advantageous  to  write  the  equations  in  a  non-dimensional 
form,  as  was  done  by  the  author  (Hoff  et  al.  1954)  in  the  absence  of 
prestresses  in  the  median  plane.  A  further  simplification  was  proposed 
by  Nachbar  in  1962  who  normalized  the  coordinates  and  displacements  for 
the  case  of  uniform  axial  prestress i  For  the  present  purposes  the  most 
convenient  form  of  the  equations  is 


F 

,xx 


(ID 


(12) 


where  the  normalized  coordinates  are  defined  as 

x=  (xVa)(2E/ocl)l/2  cp=  c^(2E/ocl)l/2  (13) 


and  the  displacements  are  given  by 

u=  (u*/a)(2E/ocl)1^2  v  =  (v*/a)(2E/crcl)1//2  w  =  (w*/a) 

(14) 

In  these  expressions  has  the  value  given  in  Eq. (l)  and  p  is  the 

ratio  of  the  critical  stress  0  of  the  present  theory  to  the  classical 

cr 

value 
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(15) 


p  =  a  /a  , 
cr'  cl 


The  stress  function  F  is  implicitly  defined  "by  the  following  expressions 

for  the  membrane  stresses  acccompanying  buckling  (that  is,  not  including 

the  prestress  0=0  )»: 

r  x  cr  ' 


0  /E  =  F 
x7  ,(f(p 


a  /E 


^ - 


-  -  F 


Finally  it  is  noted  that  subscripts  following  a  comma  indicate  differ- 
entiation  and  the  biharmonic  operator  7  is  defined  as 

A  =  QdVdx4)  +  2('dk/bx%2)  +  (bk/bcpk)\  z  (11 


Derivation  of  the  Classical  Solution 

It  is  easy  to  show  that  the  solution  of  Eqs.(ll)  and  (12)  is  the 
classical  buckling  stress  of  Eq.(l)  if  the  classical  displacement  pattern 
is  assumed.  Indeed,  with 


w  -  sin  mx  sin  no 


Equation  (12)  yields 


,  2  2x 

(m  +n  ) 


~2  sin  mx  sin  n<p 


Substitution  in  Eq. (ll)  leads  to  the  condition 


Z  +  Z"1  =  2p 


where 


,  2  22,  2 
(m  Hi  )  /m 
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If  Z  is  considered  a  continuous  variable,  the  condition  of  a  minimum 
of  p  is  that  the  derivative  of  the  left-hand  member  of  Eq.  (20)  with 
respect  to  Z  vanishes..  Thus 


-2 

1  -  Z  =  0 


(22) 


whose  solution  is  the  condition  imposed  on  the  two  reduced  wave  numbers 


2  2J~  ,  2 
[m  +n  )  /m  =  1 


(23) 


Finally  substitution  in  Eq. (20)  gives 


P  -  1 


(24) 


which  is  the  classical  solution. 


Of  course,  there  must  be  an  integral  number  of  waves  around  the 
circumference,  and  an  integral  number  of  half-waves  along  the  length  of 
the  shell.  The  first  condition  can  be  written  as 

N  =  n(2E/acl)1^2  =  integer  (25) 

If  N  is  assumed.,  n  computed,  and  m  determined  from  Eq.  (23),  in 
general,  the  number  of  half-waves  in  the  axial  direction 

M*  =  ML*/fra  (26) 

where 

M=  m(2E/acl)l/2  (27) 


will  not  be  an  integer.  However,  the  value  will  be  close  enough  to  an 
integer  in  an  engineering  approximation  whenever  the  shell  is  very  long. 
This  is  the  reason  why  Eq.(l)  has  been  accepted  generally  as  the  buckling 


formula. 
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In  a  recent  paper  Mann-Nachbar  and  Nachbar  have  shown  that  a  sta-^ 
tistical  and  probabilistic  analysis  of  the  solution  for  prescribed 
■nominal  dimensions  and  prescribed  manufacturing  tolerances  in  the 
dimensions  leads  to  most  probable  wave  numbers  and  buckling  stresses 
that  are  in  reasonable  agreement  with  experimental  data  (Mann-Nachbar 
et  al.  1965)* 

When  the  number  of  half-waves  in  the  axial  direction  is  an  integer, 
obviously 

w  =  w  =0  at  x*  =  0,L*  (28) 

2  XX 

But  from  Eqs.(l6)  and  (19)  it  follows  that 

a  =  0  =  0  and  r  4-  0  at  x*  =  0.,L*  (29) 

x  (p  xcp  T  }  ' 

This  implies  that  w  s  0  and  0^  =  0  at  x*  =  0,L*  and  thus 

v  s  0  at  x*  =  0,L*  (30) 

It  can  be  concluded  therefore  that  the  classical  solution  implies  that 

v=w  =a=v=0  at  x*  =  0,L*  (31) 

jXX  X 


This  is  obviously  not  the  only  possible  generalization  to  a  cylin¬ 
drical  shell  of  what  is  usually  assumed  to  represent  simple  supports 

(w  =  w  =  0)  for  an  ordinary  beam.  The  four  fundamental  sets  of 
,y  XX 

boundary  conditions  proposed  by  the  author  (Hoff  et  al.  1964a,  1965°) 
for  the  circular  cylindrical  shell  are: 
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Evidently,  the  SS3  condition  corresponds  to  the  classical  solution. 
Experiments  performed  in  the  ordinary  tension-compression  testing  machine 
should  be  represented  by  cases  SS2  or  SS4,  depending  upon  the  friction, 
or  its  absence,  between  testing  machine  and  test  specimen. 

The  istudy  of  the  effects  of  the  boundary  conditions  on  the  buckling 
stress  began  with  the  analysis  of  the  semi-infinite  shell.  At  the  near 
end  of  the  shell,  where  the  axial  compressive  load  was  applied,  one  of 
the  four  sets  of  boundary  conditions  given  in  Eqs. (32)  was  prescribed; 

At  the  f'lr  end,  x  °°  ,  the  displacements  and  stresses  were  required  to 
remain  bounded. 

The  solution  was  assumed  in  the  form 

w  =  e^X  sin  n<p  (33) 

Substitution  in  Eqs..  (ll).  and  (12 )  yields  eight  values  for  p  ;  these 
roots  were  first  published  by  Nachbar  in  1962.  Of  the  eight,  four  roots 
have  positive  real  parts;  these  must  be  ruled  out  for  the  semi-infinite 
shell  because  of  the  conditions  at  infinity.  The  remaining  four  provide 
us  with  four  solutions  of  the  type  shown  in  Eq.(33),  each  multiplied  by 
an  unknown  integration  constant.  The  constants  must  then  be  determined 
from  that  one  of  the  four  sets  of  boundary  conditions  given  in  Eqs. (32) 
which  is  prescribed  at  x  =  0  .  Since  the  conditions  are  homogeneous, 
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a  non-trivial  solution  of  the  four  simultaneous  linear  equations  exists 
only  if  the  determinant  formed  of  the  coefficients  of  the  unknowns 
vanishes.  Solution  of  this  buckling-  condition  yields  for  cases.  SSI 


and  SS2 


(3.4) 


The  dots  indicate  terms  that,  are  of  the  second  and  of  higher  powers  of 
(h/a)  .  They  must  be  discarded  because  in  the  derivation  of  the  differ¬ 
ential  equations  similar  terms  were  disregarded.. 


The  lowest  critical  stress  corresponds  to  N  =  2  as  the  case 
N  =  1  must  be  ruled  out  since  it  represents^  rigid -body  displacement 
rather  than  an  elastic  deformation  of  the  cross-section  of  the  shell. 
With  K  =  2  ,  Ea. (34)  becomes- 

p  -  0.5  +  1.2l(h/a)  +  .  .  .  (35) 

Since  h/a  «  1  ,  the  critical  stress  ratio  is  one-half  in  a  first 
approximation . 


Unfortunately  the  Donnell  equations  are  reliable  only  when 
»  1  ,  as  has  already  been  mentioned.  It  would  appear  therefore 
that  Eqs.(34)  and  (35)  cannot  be  trusted.  But  this  is  not  the  case, 
because  the  actual  value  of  N  is  not  important  for  very  thin-walled 
shells.  For  instance.,  when  a/h  =  1000  ,  the  assumption  of  N  =  2 
yields  p  =  0.50121  .  The  Donnell  equations  are  already  satisfactory 
when  N  =  4  .  But  with  this  value  Eq.(34)  gives  p  =  0.50484  which 
is  the  same  result  for  all  practical  purposes.  The  fact  that  the 
buckling  stress  is  insensitive  to  N  indicates  that  Eq. (35)  is  accept¬ 
able  even  though  it  was  derived  from  Donnell's  equations  for  N  =  2  . 
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This  conclusion  was  confirmed  -when,  the  author  and  his  doctoral 
student  T.  C.  Soong  recalculated  the  critical  stress  values  on  the  "basis 
of  Sanders'  differential  equations.  As  Fig.  %  shows,  the  difference 
between  the  results  based  on  Donnell's  and  Sanders'  equations  is  hardly 
noticeable  (Hoff  et  al.  1964b ^  1965b}.  On  the  other  hand,  the  Sanders 
equations  are  considered  to  be  the  best  fifst-approximatioh  theory  today, 
and  they  are  certainly  valid  when  N  =  2  (Sanders  1959) •  The  drawback 
of  Sanders'  equations  is  that  their  characteristic  roots  cannot  be  given 
in  closed  form.  Thus  Soong  had  to  evaluate  the  critical  stresses  with 
the  aid  of  a  digital  computer. 

Shells  of  Finite  Length 

On  the  other  hand,  the  author  has  succeeded  (Hoff  1964a,  1965c)  in 
obtaining  closed-form  solutions  for  the  buckling  stresses  of  shells  of 
finite  length  when  both  circular  edges  are  simply  supported  in  accord¬ 
ance  with  the  conditions  SSI.  In  this  case  all  the  eight  solutions  of 
the  type  given  in  Eq. (33)  must  be  retained,  but  the  eight-by-eight 
buckling  determinant  breaks  up  into  two  four-by-four  determinants,  one 
of  them  defining  symmetric  buckling,  and  the  other  antisymmetric  buckling. 
Of  the  two  types  of  buckling,  the  one  whose  deformations  are  symmetric 
with  respect  to  the  plane  perpendicular  to  the  axis  and  situated  half¬ 
way  between  the  two  end  planes  of  the  shell  yields  the  lower  buckling 
stress.  For  engineering  purposes  the  value  of  p  can  be  taken  as  l/2 
for  all  values  of  the  length  L*  of  the  shell  (see  Fig.  7)* 

On  the  basis  of  the  closed-form  solutions  (Hoff  et  al.  1964a,  1965c) 
and  (Hoff  1964a,  1965c)  as  well  as  the  digital  computer  solutions  (Hoff 
et  al.  1964b,  1965b)  it  can  be  stated  that  shells  of  all  lengths  whose 
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edges  are  free  to  rotate  (  w  not  restricted)  and  are  free  to  displace 

;  ’  ’  JX 

in  the  circumferential  direction  (  v  not  restricted)-  buckle  at  one- 
half  the  classical  critical  value  of  the  axial  compressive  stress, 
except  that  the  critical  stress  is  even  lower  when  the  edges  are  entirely 
free.  In  the  latter  case  p  =  0.5  for  axisymmetric  buckling,  and  0.38 
for  multi lobed  buckling  for  the  semirinf inite  shell.  For  short  shells 
p  can  be  much  smaller,  as  can  be  seen  from  the  report  by  Hoff  and  Soong 
(Hoff  et  al.  19^4b,  1965b)..  In  all  other  cases  of  boundary  conditions 
the  minimal  value  of  the  critical  stress  ratio  p  is  unity  or  greater. 

A  summary  of  this  work  was  presented  to  the  Eleventh  International 
Congress  of  Applied  Mechanics  (Hoff  1964b ).; 

A  further  digital  computer  solution  of  the  classical  equations  was 
published  recently  by  Thielemann  and  Esslinger  (Thielemann  et  al.  1964) . 
For  shells  of  finite  length  they  obtained  buckling  stresses  equal  to  the 
one  given  by  the  classical  formula  of  Eq.(l)  when  the  length  of  the  shell 
was  equal  to  or  greater  than  the  natural  wave  length  of  the  shell  and 
the  boundary  conditions  were  those  designated  in  this  paper  by  the  code 
symbols  SS3  and  SS4.  The  same  results  were  obtained  for  the  cases  RF2, 
RF3  and  RF4,  where  the  symbols  RF  indicate  rigid  end  fixation  and  the 
numerals  have  the  same  meaning  as  in  Eqs,(32).  When  the  length  of  the 
shell  was  decreased  further,  the  buckling  stress  increased.  The  same 
effect  of  the  length  was  also  found  to  exist  for  shells  reinforced  longi¬ 
tudinally  or  circumferentially. 
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Recent  Work,  by  Ohira 

The  results  quoted  are  in  complete  agreement  with,  those  obtained 
by  Ohira.  In  a  recent  private  communication  to  the  author,  Professor 
Ohira  states:  "I  got  the  original  idea  of  attempting  local  buckling 
theory  when  I  crushed  a,  number  of  Coca  Cola  paper  cups  in  a  do-it- 
yourself  laundry  during  the  time  I  was  staying  at  Purdue  in  1958"* 

His  efforts  to  develop  a.  theory  were  successful  and  in  19ol  he  presented 
his  first  paper  on  the  low  buckling  stresses  of  semi-infinite  shells 
with  modified  simple -support  conditions  (5hira-  1961).  The  following 
year  this  work  was  enlarged  (Ohira  1962),  and  in  1963  a  detailed  paper 
on  the  subject  was  presented  to  the  Fifth  International  Symposium  on 
Space  Technology  and  Science.  The  .solutions'  were  calculated  with  the 
aid  of  a  digital  computer  and  detailed  diagrams  were  presented  for  the 
SS3  and  SS.4  cases  as  well  as  for  two  cases  of  a  free  edge;  in  one  of 
these  the  u  displacement  was  kept  constant,  as  in  the  ordinary  testing 
machine,  while  in  the  second  all  displacements  were  free.  The  former 
resulted  in  p  =  0.5  ,  and  the  latter  in  p  =  0.38  .  It  is  noted  that 
unpublished  calculations  carried  out  at  Stanford  University  confirm  the 
value  0.5  for  the  first  of  the  two  free  edge  conditions.  In  all  the 
numerical  work  v  was  taken  as  0.3  and  a/h  as  300  . 

As  yet  unpublished  reports  by  Ohira,  whose  material  has  already 
been  presented  at  open  meetings  (Ohira  196h,  1965)5  deal  with  the 
buckling  of  shells  of  finite  length.  In  particular,  in  the  paper  of 
1965  the  edge  conditions  are  assumed  to  be  more  severe  at  one  end  than 
at  the  other.  The  result  is  a  satisfactory  decrease  in  the  number  of 
waves  at  buckling  and  in  the  buckling  stress  as  the  length  of  the  shell 


is  increased. 
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The  Question  of  Priority  '  . 

The  question  of  -priority  in  the  discovery  of  scientific  information 
is  always  an  interesting  one,  although  perhaps  not  an  important  one., 

It  appears  from  the  preceding  paragraphs  that  Lorenz  was  the  first  to 


calculate  the  critical  stresses  of  axially  compressed  circular  cylindrical 


shells  for  both  the  axisymmetric  and  the  general  eases  of  buckling.  Of 
course,  his  results  were  not  quite  accurate,  and  were  later  improved  upon 
by  Timoshenko  in  the  case  of  axisymmetric  buckling,  and  by  Southwell  and 
Timoshenko  in  the  general  case.  The  classical  formula  of  Eq.  (l)  was 
apparently  first  derived  by  Timoshenko. 


The  discovery  of  the  existence  of  other  solutions  of  fhe  classical 
equations  yielding  buckling  stresses  smaller  than  the  classical  one  was 
made  independently  by  Ohira  and  by  the  author  and  his  collaborators.  The 
first  oral  presentation  of  the  discovery  was  made  by  Hoff  at  a  symposium 
honoring  Dr.  von  Karman  in  Washington  on  May  11,  1961.  The  text  of  the 
talk  was  published  and  distributed  to  a  limited  number  of  recipients 
(about  200)  as  a  Stanford  University  report  in  August  1961;  the  volume 
honoring  Dr.  von  Karman,  on  the  other  hand,  appeared  only  at  the  end  of 
1962.  In  a  similar  manner,  Ohira 's  first  oral  presentation  was  made  to 
the  Eleventh  Japan  National  Congress  of  Applied  Mechanics  sometime  in 
1961  but  the  proceedings  of  the  congress  appeared  in  print  only  in  1962 
or  1963.  There  were  significant  differences  between  the  two  solutions. 
Ohira  dealt  with  the  general  type  of  buckling  of  semi-infinite  shells 
with  a  simply  supported  edge  while  the  solution  by  the  present  author 
was  valid  for  the  axisymmetric  buckling  of  semi-infinite  shells  with  a 
perfectly  free  edge.  Moreover,  Hoff's  solution  was  in  closed  form, 
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while  Ohira  had  obtained  his  results  with  the  aid  of  a  digital  computer,. 
Hoff's  result  was  generalized  in  1962  by  Nachbar  and  Hoff  to  hold  for 
the  general  case  of  buckling. 

Up  to  this  point#  Ohira  and  Hoff  had  worked  independently#  without 
knowledge  of  each  other's  work.  Subsequently  Hoff  acknowledged  that  his 
closed-form  solutions  for  the  simply  supported  edge  were  undertaken  after 
he  had  received  information  of  Ohira-' s  digital  computer  solution. 
Similarly#  in  his  paper  presented  to  the  Fifth  International  Symposium 
on  Space  Technology  and  Science  in.  Ig63#  Ohira  acknowledged  that  his 
latest  results  dealing  with  free  edges  were  obtained  after  he  had  received 
information  on  Hoff's  and  Nachbar ' s  work. 

For  a  clarification  of  these  problems  the  dates  of  first  presen¬ 
tation  and/or  publication  of  new  results  have  been  collected  in  Table  1. 

It  is  well  to  remember#  however#  that  proceedings  of  conferences  generally 
appear  ;a  year  or  more  after  the  date  of  the'  conference . 

Significance  of  New  Solutions 

It  is  equally  incorrect  to  overestimate  or  to  underestimate  the 
significance  of  the  new  solutions  of  the  small-displacement  equations. 

From  the  theoretical  standpoint  they  are  most  interesting  as  they  show 
that  solutions  other  than  the  classical  one  exist  for  the  problem  of 
the  buckling  of  axially  compressed  thin-walled  circular  cylindrical 
shells  even  though  the  problem  is  defined  by  linear  equations.  As  a 
matter  of  fact  the  only  distinction  of  the  classical  solution  is  that 
it  is  the  easiest  one  to  obtain. 
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From  a  practical  standpoint  the  hew  solutions  show  that  a  thin 
shell  placed  between  the  platens  of  a  testing  machine  can  buckle  anywhere 
between  p  =  0,5  and  1.0  ,  depending  on  the  friction  between  specimen 
and  testing  machine.  This  fact;  certainly  helps  to  explain  the  large 
scatter  observed  in  the  test  results. 

If  the  design  engineer  maintains  that  he  already  has  enough  empir¬ 
ical  information  on  the  practical  buckling  stresses  of  thin  shells  to 
make  him  uninterested  in  refinements  of  the  theory,  one  can  counter  his 
criticism  by  saying  that  the  empirical  information  is  only  on  shells 
already  built,  and  not  on  new  types  of  shells  to  be  constructed  in  the 
future.  Improvements  to  be  brought  to  shell  design  in  the  future  can 
be  evaluated  in  advance  only  if  a  complete  and  reliable  theory  lias  been 
established,  and  the  new  solutions  of  the  classical  linear  shell 
equations  make  a  contribution  to  this  goal. 

Of  course,  most  shells  in  engineering  have  their  edges  attached 
to  other  shells,  or  to  reinforcing  rings.  It  is  almost  obvious  that 
such  rings  can  provide  an  almost  perfect  restriction  of  the  circumfer¬ 
ential  displacements.  This  was  shown  to  be  the  case  in  a  recent  pub¬ 
lication  by  Almroth  (1965b)'.  On  the  other  hand,  as  yet  unpublished 
calculations  carried  out  at  Stanford  University  indicate  that  a  relax¬ 
ation  of  the  boundary  restraint  over  the  distance  of  a  single  half  wave 
can  reduce  the  buckling  stress  significantly.  Such  a  relaxation  can  be 
the  consequence  of  a  broken  rivet  or  a  poor  bond. 


THE  DISCOVERIES  OF  VON  KARMAN  AND  HIS  COIiABORATORS 

Between-  1939  and  19^1 ,  Dr.  Theodore  von  Karman  ansi  the  very  capable 
and  enthusiastic  group  of  students  and  research  men  gathered  around  him 
in  the  Guggenheim  Aeronautic  laboratory  of  the  California  institute  of 
Technology  laid  the  foundations  of  the  large -displacement  theory  of  the 
buckling  of  thin  shells  (von  Karman  et  al.  1939;  19^0,  19^-1;  Tsien  19^2a^b). 
Karman  was  puzzled  on  the  one  hand  by  the  large  difference  between  the 
complete  predictability  of  the  buckling  stress  of  rods  and  thin  plates, 
and  on  the  other  hand  by  the  very  substantial  difference  between  theoret¬ 
ical  and  experimental  values  of  the  stress  at  which  thin  shells  collapsed. 

He  also  observed  that  the  development  of  buckles  and  bulges  was  gradual 
with  elastic  rods  and  plates  while  it  was  sudden,  and  even  explosive,  with 
shells.  After  buckling,  columns  and  plates  continued  to  carry  the  buckling 
.load,  and  were  even  capable  of  supporting  further  increased  loads,  but 
with  shells  the  load  supported  after  buckling  always  dropped  to  a  fraction 
of  the  buckling  load.  Yet  the  equations  defining  the  equilibrium  and  the 
stability  of  rods,  plates  and  shells  were  all  based  on  the  same  well-proven 
hypotheses  of  the  theory  of  elasticity. 

Sudden  buckling  with  a  drop  in  the  load  carried  had  been  observed 
by  von  Karman  much  earlier  (von  Karraan  1910a)  in  the  case  of  short  columns. 
The  phenomenon  had  been  explained  by  him  completely  as  an  interaction 
between  the  inelastic  behavior  of  the  short  column  and  the  elasticity  of 
the  testing  machine  (see  also  Hoff  1961).  But  von  Karman  was  reluctant 
to  accept  a  similar  explanation  for  the  puzzling  behavior  of  the  shell, 
because  very  thin-walled  shells  appeared  to  recover  completely  after 


removal  of  the  load  and  thus  they  snowed  the  same  phenomenon  without 
plastic  deformations . 

The  clue  to  the  puzzle  seemed  to  lie  in  the  form  of  the  solution 
of  the  classical  equations,  as  given  in  Eq.(l8)„  Evidently  the  expres¬ 
sion  forthe  radial  displacement  w  remains  valid  if  the  right-hand 
member  of  the  equation  is  multiplied  by  -  1  j  it  represents  displace¬ 
ments  of  equal  magnitude  in  the  inward  and  outward  directions.  But 
real  shells  tested  in  the  laboratory  always  show  a  preference  to  buckle 
inward,  and  their  displacements  in  the  outward  direction  are  much 
smaller  than  those  oriented  inward.  It  appears  therefore  that  the 
classical  equations  fail  to  represent  properly  the  difference  between 
the  inward  and  outward  directions  in  the  case  of  shells  while  they  are 
perfectly  satisfactory  in  the  case  of  rods  and  plates.  This  difference 
must  be  a  consequence  of  the  curvature  of  the  shell  because  in  all 
other  respects  the  basic  hypotheses  upon  which  stability  theory  is  built 
are  equally  valid  for  rods  and  plates  the  one  hand,  and  shells  on 
the  other. 

The  only  manner  in  which  this  shortcc  ing  of  shell  theory  could  be 
remedied  was  to  add  terms  to  the  equations  th^t  were  non-linear  in  the 
displacements.  But  relatively  simple  equations  containing  such  terms 
had  already  been  developed  by  a  former  collabc  tor  r C  von  Karman, 

Lloyd  H.  Donnell,  who  had  worked  in  the  Guggenhei  Aeronautic  Laboratory 
of  the  California  Institute  of  Technology  betwe  n  1930  and  1933.  To 
the  linear  terms  of  his  small  displacement  equations  (Donnell  1933) 
which  were  to  become  famous  later,  Donnell  added  the  non-linear  terms 


contained  in  the  von  Karman  large-deflection  plate  equations  (von  Karman 
1910b)*  to  obtain  what  are  generally  known  today  as  the  von  Karman-Donnell 
large -displacement  shell  equations  (Donnell  1934).  The  trouble  with 
these  equations  is  of  course  that  they  cannot  be  solved  rigorously  because 
they  are  not  only  non-linear.,  but  also  of  a  very  high  order  (the  eighth). 

In  his  paper  of  1934  Donnell  showed  an  unusual  amount  of  ingenuity 
and  engineering  insight  when  he  introduced  for  the  radial  displacements 
at  buckling  expressions  equivalent  to  those  that  are  considered  today 
the  most  important  ones  in  the  representation  of  the  shape  of  the  buckles': 


w*  -  A.  cos  (,7rx*/L*)  cos  (iry*/L*)  +  7*  cos  (2.7rx*/L*)  (36) 

xx  x  y  x 


where  y*  =  a<r*  is  the  circumferential  coordinate,.  L*  and  L*  are 
T  x  y 

the  half-wave  lengths  in  the  axial  and  circumferential  directions ,  and 


and  A^q  are  constants.  After  substitution  of  this  expression  in 
his  compatibility  equation,  he  obtained  a  rigorous  solution  for  the 
stress  function.  He  used  the  principle  of  virtual  displacements  to 
calculate  the  load  necessary  to  maintain  equilibrium  at  any  value  of 
the  amplitude  of  the  buckles. 


Donnell  also  assumed  that  the  shell  was  inaccurately  manufactured. 
For  the  small  initial  deviations  from  the  exact  shape  he  introduced  an 
expression  of  the  type  represented  by  Eq.(36)  with  the  values  of  the 
constants  somewhat  complicated  functions  of  the  nominal  geometry  of  the 
shell.  The  shell  was  assumed  to  fail  when  the  stress  at  the  crest  of 


*  Apparently  Donnell  was  not  aware  of  the  existence  of  these  equations 
when  he  derived  his  large -displacement  theory. 
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the  buckles  -reached  the  yield  stress  of  the  material  of  which  the  shell 
was  constructed,  Donnell  ended  his  paper  by  proposing  a  semi -empirical 
design  formula,;  based  on  his  analysis  and  test  results,,  in  which  three 
material  constants  appeared,  namely  Young's  modulus,  Poisson's  ratio, 
and  the  yield  stress, 

The  first  innovation  in  the  analysis  by  von  Kafman  and  Tsien  (von 
Karman  et  al«  194l)  was  the  introduction  of  a  buckle  pattern  based  on 
visual  observation  of  the  shells  after  failure.  The  pattern  was  defined 
by  the  following  expression  for  the  radial  displacement: 

w*  =  A00  +  cos  (mx*/a)  cos  (ny*/a)  +  A2Q  cos  (2mx*/a) 

+  ‘.AQ2  cos  (2ny*/a)  (37) 


The  first  term  in  the  right-hand  member  represents  the  classical  small- 
displacet.cnt  solution.  The  sum  of  the  second  and  third  terms  with 
Aog  =  Ag0  defines  deformations  quite  similar  to  those  of  the  middle 
portion  of  the  test  specimen  of  Fig,  4,  but,  of  course,  the  mathematical 
expressions  imply  a  continuous  set  of  buckles  covering  the  entire  surface 
of  the  shell.  Such  a  pattern,  denoted  the  diamond  pattern,  can  be  real¬ 
ized  only  if  a  close-fitting  mandrel  is  placed  inside  the  shell  (Horton 
et  al,  1965)  (see  Fig.  8).. 

The  compatibility  equation 


(l/E)V4F  =  -(l/a)w*x*x*  +  (w*x*y*)2  "  w*x*x*w*y*y* 


was  solved  rigorously  after  the  expression  of  Eq.(37)  was  substituted 
for  w*  in  the  right -land  member.  The  second  of  the  von  Karman-Donnell 
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equation  is* 


<D/h)vV  =  -  dtf»  .  -  +  F 


,x*x*  T  r  ,y*y*w,x*x*  ”  cr  ,x*y*w,x*y* 


2F 


+  F,x*x*v,y*y*  + 


- ,x*x* 


.(39) 


where  r  a  is  the  initial  axial  compressive  stress  causing  buckling  and 
F  is  a  stress  function  from  which  the  additional  membrane  stresses 
accompanying  buckling  can  be  calculated.  This  equation  was  not  solved 
directly,  but  instead  the  direct  approach  of  the  variational  calculus 
was  used  when  the  total  potential  energy  of  the  system  was  minimized  with 
respect  to  the  three  independent  displacement  parameters  A ^  and 

A20  =  A02  * 

From  the  three  conditions  of  a  minimum  of  the  total  potential  energy 
the  three  constants  could  be  calculated,  and  thus  the  value  of  the  initial 
compressive  stress  -  a  that  corresponds  to  equilibrium  could  be  deter¬ 
mined  for  prescribed  values  of  the  parameters  p  and  q  .  Of  these  u 
was  defined-  as  the  ratio  of  the  wave  lengths 


U  = 


(40) 


and  q  was  proportional  to  the  square  of  the  number  of  waves  around  the 
circumference 

1  =  n2(h/a)  (4l) 

In  one  set  of  curves  representing  the  results  of  the  calculations 
U  was  arbitrarily  taken  as  one  because  the  wave  length  ratio  was  found 
to  be  close  to  unity  in  specimens  after  they  had  buckled  in  the  testing 
machine.  The  parameter  of  the  family  of  curves  was  q  ,  and  the  minimal 


*  This  is  in  essence  the  form  in  which  the  equation  was  given  by  Kerapner. 
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absolute  value  of  the  Stress  was 

-  °m±n  =  0.19^3 (h/a)  (42) 

The  relationship  between  the  stress  and  the  shortening  eL*  of  the 
distance  between  the  two  ends  of  the  cylindrical  shell  was  also  cal¬ 
culated  and  plotted.  The  curve  obtained  differed  little  from  the  one 
labeled  "Case  1"'  in  Fig.  9>  which  was  obtained  by  Kempr.er  in  1954. 

In  Fig-  %  the  stress  a  is  the  axial  stress,  and  it  is  considered 
positive  when  compressive.  The  abscissa  is  the  average  compressive 
strain  e  multiplied  by  the  ratio  a/h  .  The  gap  between  the  straight- 
line  and  the  curved  portions  of  the  diagram  was  not  filled  in  because 
computation  of  the  intervening  unstable  states  of  equilibrium  involves 
great  difficulties.  In  their  1941  paper,  von  Karman  and  Tsien  also 
calculated  the  equilibrium  curve  for  p  «  l/2  and  obtained  a  minimum 
for  o  which  was  negative  indicating  tension.  They  attributed  this 
unlikely  result  to  the  inaccuracies  of  their  analysis.  As  a  matter  of 
fact,  they  were  very  modest  about  their  contribution  to  science  and 
stated  that  their  rough  first  approximation  to  the  true  solution  of  the 
problem  would  have  to  be  replaced  by  a  much  more  rigorous  solution. 

Yet  von  Karman  and  Tsien  accomplished  a  great  deal.  They  discovered 
the  existence  of  three  states  of  equilibrium  corresponding  either  to  a 
prescribed  displacement  of  the  loading  head  of  the  testing  machine 
(loading  in  a  rigid  testing  machine),  or  to  a  prescribed  value  of  the 
load  (so-called  dead -weight  loading).  They  conjectured  that  in  the 
first  part  of  the  loading  process  the  states  of  stress  and  shortening 
would  follow  the  straight-line  portion  of  the  diagram  (Fig.  9)  which  is 
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stable  in  the  presence  of  infinitesimal  disturbances;  this  would  happen, 
however ,  only  if  the  shell  were  perfect  geometrically  and  if  it  had  been 
built  of  a  completely  homogeneous  and  isotropic  linearly  elastic  material* 

In  the  presence  of  initial  deviations  from  uniformity  the  maximum 
value  of  the  compressive  stress  would  be  smaller  than  that  indicated  by 
the  letter  C  in  Fig.  9>  and  the  difference  between  the  numerical  value 
0.605  for  (0/ E)(a/h)  and  the  experimental  value  would  increase  with 
increasing  values  of  the  initial  deviation.  However ,.  the  random  nature 
of  the  initial  deviations  makes  it  very  difficult  to  evaluate  the 
practical  maximal  value  of  the  buckling  stress.  Hence  von  Karman  and 
Tsien  suggested  that  for  design  purposes  the  engineer  use  the  minimal 
value  of  the  equilibrium  stress,  that  is  0.19^E(h/a)  from  an  unbuckled 
state  corresponding  to  a  somewhat  higher  value  of  the  stress  the  cylin¬ 
drical  shell  would  jump  into  a  state  of  large  displacements,  and  the 
minimal  stress  just  quoted  could  well  serve  as  a  lower  limit  to  the  stress 
at  which  the  jump  could  take  place. 

Von  Karman  and  Tsien  also  concluded  that  the  elasticity  of  the 
testing  machine  would  have  a  significant  effect  on  the  stress  at  which 
the  jump  takes  place  and  that  disturbances  of  the  test,  such  as  vibrations 
of  the  foundation  of  the  testing  machine,  would  be  an  important  contrib¬ 
uting  factor  to  the  early  failure  of  the  specimen.  These  two  conclu¬ 
sions  were  to  be  proved  incorrect  by  later  investigators. 

The  some  can  be  said  of  the  ingenious  Tsien  criterion  proposed  in 
19^2  (Tsien  1942b).  On  the  basis  of  a  detailed  and  completely  rigorous 
analysis  of  a  non-linear  model  of  a  shell,  namely  a  column  supported 
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laterally  by  an  arbitrary  number  of  nonlinear  springs  (Tsien  1942a), 

;  Tsien  suggested  that  the  minimal  equilibrium  stress  of  Fig.  9  "be  replaced 
as.  the  lower  bound  Of  the  practical  buckling  stress  by  that  particular 
stress  value  at  which  the  strain  energy  before  buckling  is  equal  to  the 
strain  energy  after  buckling,  in  a  test  in  a  perfectly  rigid  testing 
machine.  In  a  so-called  dead-weight  test,  the  total  potential  energy 
takes  the  place  of  the  strain  energy...  Tsien  realized  that  he  was.  replac¬ 
ing  a  lower  bound  by  another  lower  bound.  However,  Tsien' s  lower  bound 
was  higher,,  and  thus  closer  to  the  empirical  buckling  stress,  than  the 
earlier  one,  and  the  scanty  experimental  data  against  which  the  Tsien 
criterion  was  checked  indicated  satisfactory  agreement  between  theory 
and  experiment.  Incidentally,  a  slightly  simpler  non-linear  model  than 
the  one  studied  by  Tsien  had  been  analyzed  earlier  by  H.  L.  Cox  in  1940. 

Systems  for  which  the  Tsien  criterion  is  a  poor  approximation  were 
mentioned  by  Fung  and  Sechler  in  a  rather  complete  survey  dealing  with 
the  instability  of  shells  (Fung  et  al.  i960)  and  presented  at  the  First 
Symposium  on  Naval  Structural  Mechanics  held  at  Stanford  University  in 
1958*  Much  recent  experimental  evidence,  to  be  discussed  later,  also 
shows  that  the  proper  answer  to  the  question  of  the  practical  buckling 
stress  of  thin  shells  is  not  furnished  by  the  Tsien  criterion. 


-  32  - 


jry'spww*- 


TEDS  YOSHIMURA  BUCKLING  PATTERN 

Although  the  investigations  of  von  Karraan  and  his  collaborators 
have  resulted  in  the  discovery  of  the  physical  and  mathematical  reasons 
for  the  perplexing  behavior  of  the  axially  compressed  cylindrical  shell, 
it  was  left  to  Yoshimaru  Yoshimura,  an  imaginative  professor  in  the 
Aeronautical  Research  Institute  of  Tokyo  University,  who,  unfortunately, 
died  relatively  young,  to  find  the  geometric  reason  for  this  physical 
behavior.  Yoshimura  proved  in  a  Japanese  paper  published  in  1951  that 
the  middle  surface  of  the  circular  cylindrical  shell  is  developable  into 
a  polyhedral  surface  consisting  of  identical  plane  triangles;  such  a 
surface  is  shown  in  Fig'.  10.  His  work  was  republished  in  English  by 
the  National  Advisory  Committee  for  Aeronautics  in  1955- 

The  shell  can  therefore  be  transformed  into  such  a  polyhedral 
surface  without  stretching  its  middle  surface,  that  is  without  causing 
any  membrane  stresses  to  develop.  Small  bending  stresses  are  required, 
of  course,  to  eliminate  the  initial  curvature  of  what  are  the  plane 
triangles  after  buckling,  and  the  curvature  becomes  infinite  along  the 
edges  of  the  triangles  which  form  the  ridges  of  the  polyhedron.  It  is 
not  obvious  whether  in  the  limit  as  h/a  approaches  zero  the  work 
necessary  to  produce  the  infinite  curvature  along  the  ridges  of  a 
perfectly  elastic  shell  is  finite  or  infinite,  but  for  an  ideally 
elastic -plastic  material  certainly  a  finite  amount  of  work  suffices  to 
develop  the  ridges.  But  the  bending  stiffness  of  the  thin  wall  of  the 

3 

shell  is  proportional  to  h  while  its  extensional  stiffness  is  pro¬ 
portional  to  h  ■  Hence  a  practical  shell  is  likely  to  have  a  tendency 
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to-  avoid  extensional  deformations  more  and  more  as  its.  thickness  is 
decreased.  For  this  reason  very  thin  shells  can  "be  expected  to  buckle 
in  accordance  with  the  Yoshimura  pattern  while  thicker  ones  should  have 
more  ample  curvature  along  the  ridges. 

This  conclusion  is  borne  out  by  experiment  except  for  one  important 
modification:  the  diamond-shaped  buckles  of  Yoshimura  appear  only  in  one' 
or  two  rows  rather  than  cover  the  entire  surface  of  the  shell  as  can  be 
seen  from  Fig.  4.  This  difference  must  be  a  consequence  of  the  con¬ 
ditions  at  the  boundaries  because  the  pure  Yoshimura  pattern  is  incom¬ 
patible  with  the  circular  edge  of  the  cylindrical  shell. 

The  Yoshimura  pattern  was  discovered  independently  by  Kirste  in 
1954.  It  was  also  enthusiastically  adopted  and  studied  by  Ponsford  in 
the  Guggenheim  Aeronautic  Laboratory  of  the  California  Institute  of 
Technology  (Ponsford  1953)* 
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FURTHER  DEVELOPMENT  OF  THE  LARGE -DISPLACEMENT  THEORY 

During  the  quarter  of  a  century  that  has  passed  since  the  publi¬ 
cation  of  the  v.on  Karman-Tsien  paper  of  19^1,  most  of  the  advance  in  our 
understanding  of  the  buckling  of  circular  cylindrical  shells  subjected 
to  uniform  axial  compression  has  been  achieved  through  investigations 
using  the  von  Karman-Donnell  equations  and  the  techniques  developed  in 
that  pape-r  A  number  of  corrections  and  improvements  were  made  by 
Leggett  and  Jones  in  19^2  (but  due  to  war  conditions  “their  paper  was 
distributed  widely  Only  in  194-7),  by  Michielsen  in  1948,  and  by  Kempner 
a  doctoral  student  of  the  author,  in  1954.  In  particular,  the  total 
potential  energy  was  minimized  with  respect  to  the  parameters  q  and 
t)  defined  in  Eqs.(40)  and  (4l) .  This  minimization  showed  that  the 
buckled  state  found  by  von  Karman  and  Tsien  for  tension  (for  q  =  l/2  ) 
was  not  a  state  of  equilibrium  and'  thus  eliminated  an  inconsistency  from 
the  theory.  The  results  of  Kempner1 s  analysis  are  shown  in  Fig.  9  as 
the  curve  labeled  "Case  1". 

The  analysis  was  extended  to  orthotropic  shells  in  a  paper  presented 
by  Thielemann  at  the  Durand  Centennial  Conference  (Thielemann  i960)  and 
the  results  of  the  calculations  were  compared  with  experiments  carried 
out  by  Thielemann  with  extreme  care.  This  work  was  continued  by 
Thielemann  in  a  report  to  a  NASA  conference  held  at  Langley  Field  in 
1962  in  which  he  objected  to  the  minimization  of  the  total  potential 
energy  with  respect  to  the  wave  lengths  because  evidently  only  integral 
numbers  of  waves  can  occur  in  the  shell.  It  is  not  clear,  however,  from 
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the  very  concise  paper  what  procedure  he  intro to  replete-  this 
minimization.  ' 

Thielemann  used  electronic  digital  comouters  to  solve,  ir.  an 
approximate  -manner*  the  von  Knrroan -Donnell  equations.  The  use  ot  tne 
digital  computer  was  exploited  even  more  completely  by.Aimroth  who 
investigated  many  combinations  of  the  various  terms  in  the  series 


00  CO 

=  ^  ^  cos  (j7Tx*/L*)  cos  (kvao^/L*) 


j-0  k=Q 


in  order  to  obtain  the  minimum  of  the  total  potential  energy  The 
curves  labeled  "Case  2"  and  "Case  3"  in  Fig.  9  represent  Airmotn's 
results;  in  the  calculations  of  the  former  only  the  coefficients  A^  , 

A11  3  A22  5  A20  an<^  \o  were  assuine<*  t)e  different  from  zero 

while  in  the  latter  the  coefficients  A_,  and  A,^  were  also  included 

33  60 

In  order  to  make  the  computational  work  tolerable,  all  the  coefficients 
not  listed  were  assumed  to  be  zero. 


It  can  be  seen  that  for  a  fixed  value  of  ea/h  the  equilibrium 
stress  decreases  as  the  number  of  terms  considered  is  increased.  The 
important  question  to  ask  is  therefore  where  the  limiting  curve  is 
situated  when  the  number  of  terms  considered  is  further  increased  and 
made  to  approach  infinity.  Almroth  felt  that  his  nine-term*  approxi¬ 
mation  (not  shown  in  the  figure)  approached  closely  enough  the  limiting 
curve;  he  was  unable  to  obtain  significant  changes  by  selecting  different 
coefficients,  or  considering  additional  ones.  Moreover,  his  results 
were  in  excellent  agreement  with  those  obtained  both  theoretically  and 
experimentally  by  Thielemann  (1962). 
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The  author  of  the  present  paper  was  not  convinced  that.  Almroth's 
nine -ter#  approximation  was  a  sufficiently  good  approximation  to  the. 
limiting  curve  but  he  realized  that  it  would  be  extremely  difficult  to 
continue  the  process  by  adding  more  and  more  terms  to  the  displacement 
expression.  In  particular ,  the  replacement  of  the  products  and  powers 
of  trigonometric  terms  by  trigonometric  terms  of  multiple  angles  was 
such  a  lengthy  job  in  the  analysis  of  the  problem  that  it  was  almost 
impossible  to  avoid  errors.  For  this  reason  the  author  suggested  that 
this  work  should  be  programmed  for  the  computer;  the  method  developed 
for  this  purpose  has  been  described  in  two  reports  (Madsen  et  al.  1965b; 
Bushnell  et  al..  1965)* 

With  the  aid  of  this  computer  program  the  author  and  his  collabora¬ 
tors  (Hoff  et  al.  1965a)  succeeded  in  calculating  equilibrium  curves  on 
the  basis  of  up  to  l4-term  approximations .  Hn  manner  in.  which  the 
equilibrium  stress,  decreases  with  the  number  of  terms  considered  for  a 
fixed  value  of  ea/h  can  best  be  seen  from  the  entries  under  cases  4,  5 
and  6  in  Table  2.  The  normalized  stress  values  are  not  minimal  values 
in  this  instance,  but  they  correspond  to  ea/h  =  3.4  .  They  are  O.O856, 
0.0706  and  0.0528  when  the  number  of  terms  retained  is  8,  10  and  12. 

When  two  more  terms  were  added  to  the  series  the  value  of  (o/E)(a/h) 
became  0.0427,  which  is  the  minimum  of  the  curve  labeled  "Case  7"  in 
Fig.  9. 

*  In  the  discussion  that  follows  the  A^q  term  is  not  counted  when  the 
total  number  of  terms  is  indicated.  The  value  of  Aqq  is  obtained 
from  considerations  of  the  continuity  of  deformations,  and  not  from  a 
minimization. 
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This  value  is  substantially  smaller  than  Almrbth 1 s  0.0652;  It  Is  also 
smaller  than  the  value  of  0. 65 18  obtained,  by  Sobey  ( l$6Vb )  with  a  23-term 
approximation.  The  inclusion  of  such  a  large  number  of  terms  was  possible 
only  because  of  the  availability  of  a.  superior  computer  program.  The 
reason  why  Sobey' s  stress  value  is  higher  for  23  terms  than  the  author's 
value  for  l4  terms  is  thav.  many  of  the  coefficients  of  the  terms  retained  . 
by  Sobey  have  very  small  numerical  values.  Hence  these  terms  are  unimpor¬ 
tant  in  the  definition  of  the  displacement  pattern.  Incidentally,  Sobey 's 
paper  had  a  very  limited  distribution  and  was  unknown  to  the  author  at  the 
time  he  wrote  his  paper  jointly  with  Madsen  and  Mayers. 

The  most  interesting  result  of  the  paper  by  Hoff,  Madsen  and  Mayers 
(1965a)  is  the  observation  that  with  increasing  numbers  of  terms  retained 
in  the  expression  for  w.  the  coefficients  of  the  terms  of  the  double 
Fourier  series  approach  the  values  characterizing  the  Fourier  expansion  of 
the  Yoshimura  buckle  pattern.  At  the  same  time,  4  ,  r\  and  a  approach 

zero. 

It  appears  therefore  that  the  shell  buckles  into  an  exact  Yoshimure 
pattern,  with  a  finite  wave  length  in  the  axial  direction  but  a  vanishing 
wave  length  in  the  circumferential  direction  (p  =  0).*  At  the  same  time 
t|  =  N^h/a  approaches  zeroj  since  the  number  of  waves  around  the  circumfer¬ 
ence  cannot  be  less  than  two,  obviously  h/a  must  approach,  zero.  In  other 
words,  the  limiting  curve  obtained  when  the  number  of  terms  is  increased  and 
made  to  approach  infinity  is  a  rigorous,  but  trivial,  solution  because  it  is 
valid  only  for  an  infinitely  thin  shell.  Evidently  the  stress  under  which 
an  infinitely  thin  shell  can  be  in  equilibrium  after  buckling  is  infinitely 
small. 


* 

Another  possibility  is  a  finite  wave  length  in  the  circumferential 
direction  and  an  infinite  wave  length  in  the  axial  direction. 


* 


„ . . . 


The  solution  of  this  puzzle  was  presented  in  a  follow-up  report 
by  Madsen  and  Hoff  (1965a).*  For  a-  given  cylindrical  shell,  that  is  for 
a  prescribed  value  of  a/h  ,,  t|  cannot  assume  a  value  smaller  than 
4a/h  .  lienee  minimization  with  respect  to  rj  means  a  differentiation 
of  the  total  potential  energy  with  respect  to  tj  ,  setting  the  resulting 
expression  equal  to  zero,  and  solving  for  ,  provided  that  a  value 
equal  to  or  greater  than  4a/h  is  obtained  by  this  procedure  (and  pro¬ 
vided  „hat  the  inaccuracy  connected  with  the  replacement  of  the  integral 
values  of  N  with  a  continuous  function  is  considered  admissible).  If 
the  value  obtained  for  t)  is  less  than  4a/h  ,  it  has  to  be  replaced 
by  4a /h  .  In  this  manner  a  lower  bound  exists  for  rt  and  the  minimal 
value  of  the  postbuckling  stress  is  greater  than  zero. 

In  another  extension  of  the  large -displacement  investigations 
initiated  by  von  Karman  and  Tsien  in  1941,  the  behavior  of  initially 
slightly  inaccurate  circular  cylindrical  shells  was  studied  (see  Fig.ll). 
The  importance  of  initial  deviations  from  the  exact  shape  had  been 
recognized  much  earlier  as  it  had  already  been  studied  by  Fliigge  in 
1932  and  by  Donnell  in  1934.  But  in  1950 ,  Donnell  and  Wan  greatly 
altered  the  procedure  followed  by  Donnell  sixteen  years  earlier  and 
developed  a  new  method  of  calculation  which  was  to  be  copied  by  several 
other  investigators.  Through  a  rather  complex  reasoning,  and  on  the 
basis  of  his  broad  experience  in  engineering,  Donnell  came  to  the  con¬ 
clusion  that  the  most  dangerous  initial  deviations  of  the  middle  surface 
of  a  circular  cylindrical  shell  from  the  exact  shape  could  be  represented 
by  the  equation 
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wj/h  =  ( u/rr2 )  ( L* 1  ‘  5L*° ;  5/h2 )  f(x*,y*)  =  (Ue2/n1,5lfti2)  f(x*,y*)  (44) 


where  U  is  the  unevenness  parameter  and  f(x*,y*)  the  function  given 

in  Eq.(37).  The  additional  radial  displacements  caused  by  the  load  were 

represented  in  the  same  form  and  were  multiplied  by  an  amplification 

factor.  The  compatibility  equation  (38)  was  solved  rigorously  for  the 

stress  function  P  .  The  expressions  for  w*  .  and  F  were  then  sub- 

tot 

stituted  in  the  expressions  for  the  total  potential  energy  and  the 
expression  so  obtained  was  minimized  with  respect  to  the  amplification 
factor  and  A^q  ,  A ,  m  and  n  . 


This  implies  that  the  shape  of  the  displacements  caused  by  the 
loads  was  taken  to  be  the  same  as  the  shape  of  the  initial  deviations. 
This  is  obviously  a  restriction  on  the  generality  of  the  solution,  but 
in  view  of  the  difficulties  inherent  in  any  solution  of  the  governing 
equations  it  is  a  justifiable  one.  If  the  minimization  had  been  carried 
out  only  with  respect  to  the  amplification  factor,  the  result  could  be 
accepted  os  a  usable  approximation.  Unfortunately,  the  total  potential 
energy  was  also  minimized  with  respect  to  the  parameters  defining  the 
shape  of  initial  deviations.  This  means  that  the  system  whose  total 
potential  energy  was  minimized  was  not  defined  at  all,  but  changed  its 
initial  shape  during  minimization.  A  correct  and  complete  analysis 
rhoul';  define  the  initial  shape  by  ;  ns  of  one  coefficients  A?,  , 


‘11 
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,  A_0  and  the  additional  displacements  by  means  of  the  coefficients 


A1  -  A1  A1 
l1  ’  20  '  02 


The  minimization  should  then  be  carried  out  with 


respect  to  A^  ,  and  A^  >  and  not  with  respect  to  the  coefficients 

.o  .o  ,  .o 
An  ,  A20  and  . 
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The  .same  error  can  he  found  in  a  number  of  publications  based  on 
the  paper  by  Donnell  and  Wan;  they  are  the  articles  by  Loo  (195*0  >  Lee 
(1962)  and  Sobey  (1964b). 

The  error  was  avoided  by  Madsen  and  Hoff  who  used  a  two-term 
expression  to  define  the  shape  of  initial  deviations  and  a  three -term 
expression  for  the  additional  displacements  (Madsen  et  al.  1965a)-.  The 
minimization  was  carried  out  with  respect  to  the  three  coefficients  of 
the  additional  displacements  and  the  results  are  shown  in  Fig.  11. 

It  is  evident  from  this  figure  that  small  initial  deviations  from 
the  exact  cylindrical  shape  have  a  large  effect  upon  the  maximum  load 
carried  by  the  compressed  shell;  and  it  is  worth  noting  that  this  maximum 
load  is  the  only  quantity  that  can  be  observed  directly  in  a  compression 
test.  For  instance,,  an  initial  amplitude  of  the  nonsymmetric  deviations 
amounting  to.  one -tenth  of  the  wall  thickness  coupled  with  an  amplitude 
of  the  axisymmetric  deviations  amounting  to  one-fortieth  of  the  wall 
thickness  reduces  the  maximum  load  to  60  percent  of  the  classical  value 
calculated  for  the  perfect  shell. 

A  more  complete  calculation  by  Almroth  (1965b,  1966)  resulted  in 
somewhat  lower  maximal  values  of  the  stress. 

All  the  solutions  of  the  large-displacement  equations  quoted  assume 
that  the  shell  is  very  long  and  that  its  surface  is  completely  covered 
with  uniform  bulges  after  it  has  buckled.  Yet  Fig.  4  clearly  shows  that 
in  the  laboratory  shells  buckle  only  over  a  small  area  and  that  the 
remainder  of  their  surface  remains  smooth.  This  fact  was  already  mentioned 
by  Yoshimura  (1951)*  The  localized  buckle  pattern  was  introduced  into  an 

-  4l  - 


energy  solution  of  the  large -displacement  equations  by  Uemura,  a  former 
collaborator  of  Yoshimura.,  while  he  was  working  on  a  research  project 
at  Stanford  University.  His  solution  (Uemura  1963..,  196*0  indicates  a 
tendency  on  the  part  of  the  shell  to  prefer  local  buckles  to  uniformly 
distributed  ones,  but  the  results  are  not  really  conclusive  because  of 
the  comparatively  small  number  of  terms  retained  in  the  series  repre¬ 
senting  the  radial  displacements. 

This  chapter  would  be  incomplete  without  mention  of  the  efforts 
made  to  check  whether  the  von  Karman-Donnell  equations  are  sufficiently 
accurate  for  an  analysis  of  the  postbuckling  behavior  of  thin-walled 
circular  cylindrical  shells.  On  the  one  hand  it  is  easy  to  show  that 
the  Donnell  expressions  for  the  membrane  strain  are  completely  inadequate 
to  represent  the  inextensional  deformations  of  the  Yoshimura  pattern 
when  there  are  5  to  10  triangles  around  the  circumference  of  the  shell 
(Hoff  et  al.  1965a),  and  on  the  other  the  curvature  expressions  become 
inaccurate  and  the  stz’esses  can  exceed  the  yield  stress  of  the  material 
when  the  computations  are  carried  out  with  the  retention  of  more  and 
more  terms  of  the  infinite  series  for  a  prescribed  value  of  the  a/h 
ratio.  The  latter  two  observations  were  made  by  Mayers  and  Rehfield 
in  a  report  published  in  196*1. 

Moreover  the  tremendous  effort  made  by  many  investigators  in  the 

last  25  years  has  resulted  only  in  a  reduction  of  the  value  of  the 

coefficient  k  in  the  buckling  stress  formula  0  =  kE(h/a)  ,  but 

cr 

the  value  of  k  has  remained  a  constant,  independent  of  the  a/h  ratio. 
Experiments  show,  however,  that  k  can  be  as  low  as  0.3  when  a/h  is 
100,  and  0.0b  when  a/h  is  3000. 
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It  vas  desirable  therefore  to  investigate  the  effect  upon  k  of 
the  use  of  equations  more  accurate  than  the  von  Karman -Donnell  equations. 

This  was  done  by  Mayers  and  Rehf ield  in  the  paper  cited;  they  found, 
however,  that  the  dependence  of  k  on  the  a/h  ratio  is  negligibly 
small.  The  same  conclusion  was  drawn  by  Tsao  (1965)  snd  by  Madsen  and 
Hoff  (1965).  Unfortunately,  the  calculations  of  the  former  were  shown 
to  be  unreliable  by  Mayers  and  Rehf ield.  In  the  Madsen-Eoff  article  perfectly 
rigorous  membrane  strain  expressions  and  almost  perfectly  rigorous  cur¬ 
vature  expressions  were  developed  for  arbitrarily  large  displacements  and 
for  strains  that  are  small  compared  to  unity.  The  calculations  involved 
the  minimization  of  a  total  potential  energy  expression  containing  more 
than  12,000  terms. 

The  perplexing  conclusion  must  be  drawn  therefore  that  even  though 
displacement  patterns  can  easily  be  devised  for  which  the  Donnell  strain- 
displacement  and  curvature -displacement  relations  are  grossly  inadequate, 
and  although  these  relations  form  the  basis  of  the  von  Karman-Donnell 
large -displacement  equations,  replacement  of  these  equations  by  more 
accurate  ones  does  not  change  noticeably  the  equilibrium  states  obtain¬ 
able  from  the  equations.  The  explanation  of  the  paradox  is  probably 
that  the  procedures  used  to  solve  the  equations  always  lead  to  displace¬ 
ment  patterns  involving  so  many  waves  around  the  circumference  that  the 
shell  can  be  considered  a  shallow  one. 
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THE  EFFECT  OF  PREBUCKLING  REFORMATIONS 

It  has  already  "been  mentioned  that  in  his  fundamental  paper  of 
1932  on  the  buckling  of  cylindrical  shells,  Flugge  investigated  the 
effect  on  the  buckling  load  of  deformations  that  occur  during  the  loading 
of  the  shell  before  it  buckles.  These  deformations  arise  because  of  the 
tendency  of  the  compressed  shell  to  expand  uniformly,  and  because  of  the 
restriction  of  this  expansion  at  the  supports.  Flugge 's  study  indicated 
that  because  of  the  prebuckling  deformations  the  yield  stress,  of  the 
material  is  reached  in  the  shell  slightly  before  the  critical  value  of 
the  stress  is  reached  even  though  the  critical  stress  is  well  within  the 
elastic  limit  of  the  material. 

Recently  the  problem  was  attacked  again  by  two  investigators  who 
worked  almost  simultaneously  and  entirely  independently,  without  knowledge 
of  each  other's  work*  But  in  this  new  research  the  full  power  of  the 
electronic  digital  computer  was  used  to  obtain  the  results.  In  the 
United  States  the  veteran  investigator  Manuel  Stein  (1962,.  1964)  and  in 
Germany  the  young  research  man  G.  Fischer  (1963)  solved  first  the  rela¬ 
tively  simple  axisymmetric  problem  of  the  prebuckling  deformations.  Next 
extensive  computer  programs  were  developed  for  the  solution  of  the  lin¬ 
earized  stability  problem  of  the  deformed  and  still  axisymmetric,  but  no 
longer  cylindrical  shell. 

When  the  results  were  finally  compared,  surprisingly  Fischer 's 
buckling  stress  was  found  to  be  about  twice  as  high  as  that  obtained  by 
Stein.  The  former  showed  buckling  at  about  82  percent  of  the  classical 


critical  stress  while  the  latter's  shells  buckled  at  stresses  amounting 
to  42  to  48  percent  of  the  classical  critical  value.  The  discrepancy  was 
explained  when  the  new  solutions  of  the  classical  linear  equations  obtained 
byOhira  and  Hoff,,  described  at  the  beginning  of  this  paper,,  became  known, 
in  his  analysis  Fischer  used  the  boundary  conditions  designated  by  the 
code  symbol  SS3,  and  Stein  those  denoted  SS2  (see  Eqs.(32)). 

A,  final  comparison  of  the  two  solutions  was  made  by  Almroth  (19.65b) 
who  studied  eight  sets  of  boundary  conditions,  namely  those  indicated  by 
the  present  author  by  the  symbols  SSI  to  SS,4,  and  RF1  to  RF4.  He  con¬ 
firmed  Stein's  and  Fischer's  solutions  and  concluded  that  the  effect  upon 
the  buckling  stress  of  the  boundary  conditions  was  large,  and  that  of  the 
prebuckling  deformations  was  small; 


THE  KOITER  THEORY 

Probably  the  easiest  way  to  acquire  an  understanding  of  the  funda¬ 
mental  idea  of  the  Koiter  theory  is  to  work  out  an  example  in  some  detail, 
and  for  the  example  the  non-linear  model  of  .a  shell  analyzed  recently  by 
the  author  (Hoff  1965a)  may  well  be  chosen.  The  model  (Fig-  12)  consists 
of  two  pih-^ointed  bars  whose  common,  end  point  is  supported  laterally  by 
a  non-linear  spring,  and  whose  far  ends  are  under  the  action  of  equal 
and  opposite  forces  P  .  In  the  original  paper  the  bars  were  elastic 
and  two  linear  springs  attached  to  their  far  ends  represented  the  elas¬ 
ticity  of  the  testing  machine.  Since  the  effect  of  these  features  of 
the  model  on  the  buckling  phenomenon  is  small,  they  are  omitted  from  the 
present  analysis  in  order  to  save  space  and  effort. 

The  spring  force  S  is  characterized  by  the  equation 

S  =  2Kf (£)  (45) 

where  the  non-dimensional  displacement  £  is  defined  as 

{  «  q/h  =  (y/h)  -  (e/h)  (46) 

and  e  and  h  are  the  eccentricity  of  the  system  and  the  initial 
vertical  component  of  the  length  of  each  bar.  In  the  original  publication 
the  spring  characteristic  was  defined  as 

S  =  100q3  -  200n2  +  10 5q  (47) 

A  graph  of  this  relationship  is  shown  in  Fig.  13.  In  the  present  paper 
the  quantities  appearing  in  Eq.(45)  are  defined  as 
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Substitution  in  Eq.(  51)  yields 


X  =  P/K  (55) 

First  the  system  without  eccentricity,  the  so-called  perfect  system, 
will  be  examined.  For  such  a  system 
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L  JJ 

since 


e  =  0  =  5 


(56a) 


In  agreement  with  the  principle  of  virtual  displacements  the  first 
variation  of  the  total  potential  energy  must  vanish  for  equilibrium. 

The  variation  must  be  carried  out  with  respect  to  the  only  independent 
displacement  quantity  £  .  One  obtains 


SU  =  (dU/dC)6£  =  2hK 


f(t) 


0 


(57) 


This  equation  has  two  kinds  of  solutions.  Since  f(0)  =  0  ,  evidently 
one  solution  is 


t  =  0  (58) 

This  means  that  the  system  is  in  equilibrium  in  its  fundamental  state, 
the  initial  straight-line  configuration,  whatever  the  value  of  the  load 
factor  X  . 
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In  addition,  equilibrium  is  possible  for  certain  combinations  of 
load  and  displacement  characterized  by 

P  1/2 

x®  (iA)(i-r)  f(c)  when  :  /o  (59) 

This  equation  defines  a  buckled  state  adjacent  to  the  fundamental  state 
which  will  be  called  the  adjacent  state-. 

The  stability  of  the  fundamental  state  depends  upon  the  sign  of  the 

second  derivative  of  the  total  potential  energy..  From  Eq. (57)  ohe 

obtains : 


(afyat2)  =  2hK|T'(t)-Ui-^) 


2-3/2- 


But  for  the  fundamental  state  t;  =  0  ;  hence 


(d2u/dt2)  c=0  =  2hK[f  •  (0)  -  x] 

From  Eq.  (48)'  evidently 
f'(0)  =  1 


(61) 


(62) 


Consequently  the  fundamental  state  is  stable  when  X,  <  1  and  unstable 
when  K  >  1  .  The  critical  point  is  characterized  by 

X  =1  that  is  P  =  525  lb  (63) 

cr  cr  ' 


The  critical  point  is  a  bifurcation  point,  or  branching  point,  where 
equilibrium  configurations  adjacent  to  the  fundamental  configuration 
appear  for  the  first  time  in  the  loading  process.  The  branching  point 
is  labeled  Q  in  Fig.  14. 

Koiter  was  the  first  to  call  attention  to  the  importance  of  the 
stability  of  the  system  in  the  branching  point  itself.  There  the  second 
derivative  of  the  total  potential  energy  is  zero  and  thus  stability 


depends  on'  the  derivatives  of  higher  order.  'From  Eq.(60)  the  third 
derivative  is  easily  obtained: 


(d  u/dt°)  =  2hK 


f"(0-3 


-5/2: 


m 


Since  in  the  branching  point  £  =  0  and  X  =  1  ,  the  expression  becomes 


(a5u/ac3)t,0  =  2hKf"(o) 

X=1 


(65) 


From  Eq. (48)  one  calculates 


f  "(0)  =  2ag  =  r  4000/105  «  -  38 


(66) 


Thus 

(d3uM3)^=0  *  -  76hK  4  0  (67) 

X=1 

When  the  second  derivative  of  the  total  potential  energy  vanishes 
and  at  the  same  time  the  third  derivative  is  not  zero,,  the  system  is 
unstable  (see.,  for  instance,  Hoff  1956).  Evidently  in  such  a  case  the 
equilibrium  corresponds  to  a  rainimax,  and  the  total  potential  energy 
decreases  during  a  small  positive  excursion  if  it  increases  during  a 
small  negative  excursion,  and  vice  versa.  Koiter  has  shown  that  under 
such  conditions  the  system  is  very  sensitive  to  small  initial  deviations 
from  the  perfect  shape. 

This  sensitivity  can  be  checked  if  the  imperfect  system  is 
investigated.  The  total  potential  energy  is  given  by  Eq.(54);  its 
first  derivative  is 
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-l/2Y 

=  2hK(f  (5)  -  xf£+8)[l  +  5s  -  ( C+S)2]  >  =  0 


(63) 


With,  a  non-vanishing  eccentricity  e  =  h&  the  above  expression  of  the 
principle  of  virtual  displacements  has  only  one  solution: 


\=  a+srVoCi+t2-  (C+h)2^ 


(69) 


The  second  derivative  of  U  is 


2 stl  -2  /b  ^..2l'3^2 


d  U/df  =  2hKvf  >  ( C )  -  x(l+&2)  [l  +  o  -  (£,+&) 


J 


(TO) 


For  sufficiently  small  absolute  values  of  &  and  X  this  expression 
is  certainly  positive;  hence  the  load-displacement  curve  defined  by 
Eq. (69)  is  stable  when  the  load  is  small.  The  stability  vanishes  when 
the  second  derivative  becomes  zero.  Substitution  of  the  expression  for 
X  from  Eq. (69)  and  equation  to  zero  of  the  second  derivative  result  in 

-1 


f*(o  -(i+52)u+&)“1[i +62-  a+sr].  f(o  =  0 

When  1  and  5  are  sufficiently  small,  this  simplifies  to 


(TD 


f'(t)  -  (t+5)_1f(t)  =  0 


(72) 


In  view  of  the  graph  of  S  =  2Kf(£)  shown  in  Fig.  13  this  equation 
has  no  real  solution  when  6  and  £  are  negative.  The  curves  repre¬ 
senting  the  displacements  of  a  system  whose  eccentricity  is  negative 
are  stable  everywhere.  A  critical  point  can  exist,  however,  when  & 
and  t,  are  positive,  but  this  critical  point  is  not  a  branching  point 
but  a  limit  point,  that  is  a  maximum  of  the  load -displacement  curve. 

The  curves  shown  in  Fig.  l4  indicate  that  the  maximum  of  the  load 
reached  by  a  slightly  imperfect  system  can  be  much  lower  than  the 


-  51  - 


classical  critical  load  of  the  perfect  system.  This:  is  always  the  case 
when  the  branching  point  of  the  perfect  system  .(point  Q  in  Fig.  it)  is 
unstable; 

The  opposite  is  true  when  the  branching  point  of  the  perfect  system 
is  stable,  as  will  now  be  shown.  Let  us  attach  a  second  spring  to  the 
.joint  of  the  system,  shown  in  Fig.  12,  but  in  the  opposite  direction. 

The  horizontal  force  S'  provided  by  the  second  spring  will  then  be 

S'  =  100t]3  +  200q2  +  105n  (73) 

To  maintain  unchanged  the  classical  critical  load>  the  dimensions 
of  the  springs  will  be  reduced  until  each  provides  only  one-naif  the 
force  it  did  before.  The  combination  of  the  two  springs  will  now  be 
characterized  by 

S"  =  (1/2)  (S+S')  =  100  r] 3  +  105  q  (74) 

and  in  Eqs.(48)  the  only  change  to  be  made  is  to  write 

f(C)  =  a3C3  +  C  (75) 

The  total  potential  energy  expression  of  Eq.(56)  anu  the  expressions 
for  the  derivatives  given  in  Eqs,(57)>  (60)  and  (64)  remain  unchanged. 
Again,  the  bifurcation  of  the  equilibrium  states  occurs  at  P  =  525  lb 
and  the  fundamental  state  is  stable  below,  and  unstable  above  this  value. 
But  the  second  derivative  of  f (£,)  is  different: 

f"(0  =  6C  (76) 

In  the  fundamental  state  this  obviously  vanishes.  Hence  in  the  branching 
point 
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But  from  Eqs.(48)  and  (75), 

fm  (0  =  6a-3  =  >600(100/105)  (80) 

It  can  be  concluded  therefore  that  in  the  branching  point  the  fourth 
derivative -of  the  total  potential  energy  is.  positive,  and  thus  the 
equilibrium  of  the  branching  point  is  stable. 


The  equilibrium  states  were  also  investigated  in  the  presence  of 
small  initial  eccentricities  and  the  curves  representing  the  behavior 
of  the  system  are  shown  in  Fig.  15.  It  can  be  seen  from,  the  figure 
that  imperfections  have  no  significant  effect  upon  the  load  the  system 
can  carry. 


Figure  15  is  representative  of  the  behavior  of  a  flat  rectangular 
plate  compressed  in  its<  plane  with  its  edges  simply  supported..  After 
buckling  the  load  can  be  increased  further  and  small  deviations  from 
flatness  have  little  effect  on  the  load-carrying  capacity  of  the  plate. 
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Figure  14  is  characteristic  of  the  behavior  of'  an  axially  compressed 
thin-walled  circular  cylindrical  shell.  The  maximum  load,  it  can  carry 
is  greatly  affected  by  small  deviations  from  the  exact  cylindrical  shape 
and  the  load  drops  suddenly  when  the  critical  value  of  the  imperfect 
system  is  reached',  in  the  testing,  machine . 

The  connection  between  postbuckling  behavior  and  the  stability  of 
the  system  was  explored  in  detail  for  such  complex  systems  as  shells, 
and  criteria  for  determining1  the  stability  of  the  branching  point  were 
established  rigorously  in  a  doctoral  dissertation  written  by  Koiter  in 
1945 •  Unfortunately,  the  dissertation  was  published  -in  the  Dutch 
language  and  for  a  long  time  it  did:  not  receive  the  attention  it  merited. 
A  concise  presentation  of  the  principles  involved  was  made  by  Koiter  at 
the  Symposium  on  Non-Linear  Problems  in  Madison,  Wisconsin,  in  1963  and 
the  printing,  of  the  paper  in  the  Proceedings  of  the  -symposium  has  con¬ 
tributed  greatly  to- the  recognition  of  the  importance  of  the  theory  in 
the  analysis  of  structural  stability.  A  third  publication  by  the  same 
author  (Koiter  1963Vv)  contains  a  rigorous  solution  for  the  imperfect 
circular  cylindrical  shell. 


It  follows  from  Koiter 's  .general  theory  that  the  ratio  p  of  the 
maximum  stress  of  the  imperfect  shell  to  the  classical  critical  .stress 
of  the  perfect  shell  is  given  by  the  equation 


Pt  =  [4/27(l-V2)]  ^  (1-p) 


*(8l). 


if  the  imperfections  are  axially  symmetric  and  \|r  is  the  ratio  of  the 
amplitude  of  the  sinusoidal  initial  deviations  from  the  exact  circular 
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cylindrical  shape  and  the  wail  thickness.  Introduction  of  the  new 
symbols 


z  =  1  -  p 


C  * 


-v; 


\Jr  »  2.48i|r 


permits  the  writing  of  the  solution  of  Eq.;(8l)  in  the  form 


z  - 


a  ,  -  a 
cl  .,-cr 

a  , 
cl 


=  -  (c/2)  +C1/2  Q:  +  ( 


.-,1/2 


When  C/4  <<  1  ,  this  is  equivalent  to 
z  =  C1/2  -  (i/2)c  +  (i/8)C3/2  +  .  . 


(82) 


(83) 


(84) 


In  an  earlier  paper  (Madsen  et  al.  1965a).  it  was  proposed  that  in 
a  first  approximation  the  inital  deviation  amplitude  should  be  assumed 
to  be  proportional  to  the  radius  of  the.  shell.  Since  \|r  Is  this 
amplitude  .divided  by  the  wall  thickness.  h.  one  can.  write 


t  =  K*(a/h) 


(85) 


The  formulas  given  lead  to  reasonable'  agreement  with  experimental  data, 

as  was  indicated  by  the  author  in  his  lecture  at  the  Seventh  International 

Aeronautics  Congress  in' Paris  (Hoff  1965b),  if  the  value  of  K*  is 
,-4 


chosen  as  10 


If  one  wants  to  obtain  a  formula  valid  for  less  care¬ 


fully  manufactured  specimens  he  may  choose 


-4 


K*  =  4  x  10 


Substitutions  yield 


C  10"3(a/h) 


(86) 


(87) 
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From  Eqs. (82)  and  (83)  the  values,  of  p  become  about  0.64,  0.47,  0.38 
and  0.21  when  the  a/h  ratio  is  200,  600,  1000  and  3000. 

The  KOiter  theory  has  recently  been  taken  up  by  investigators  at 
University  College  in  London  and  at  Harvard  University  in  Boston  and 
a  nuniber  of  interesting  results  were  obtained  in  the  former  place  by 
J.  M.  T.  Thompson  (Thompson  1961,  1963,  1964)*  and  in  the  latter  by 
J,  Hutchinson  (Hutchinson  .1965) .. 


*  Part  of  Thompson's  work  was  carried  out  at  Stanford  University. 


EXPERIMENTAL  VERIFICATION 

The  development  of  the  theory  has  always  gone  hand  in  hand  with 
increasingly  careful  experimentation.  No  details  will  be  given  here  of 
experiments  conducted  to  determine  the  buckling  stresses  of  circular 
cylindrical  shells.  It  may  be  mentioned,  however.,  that  the  derivation 
of  the  classical  buckling  stress  formula  was  preceded  by  Lilly’s  experi¬ 
ments  in  1908,  and'  that  the  revival  of  interest  in  buckling  theories  in. 
the  thirties  was  paralleled  by  the  experimental  work  of  Robertson  (1928, 
1929),  Fliigge  (1932),  Wilson  and  Newmark  (1933):,  Lundquist  (1933), 
Donnell  (1934)  and  Kanemitsu  and  Nojima  (1939). 

In  the  more  recent  past  large-scale  experiments  were  carried  out 
with  specimens  of  large  a/h  ratios  by  Harris,  Suer,  Skene  and  Benjamin 
in  1957  and  Ly  Weingarten,  Morgan  and  Seide  in  1965*  Thielemann  (i960, 
1962)  also  made  a  large  number  of  tests  at  the  time  when  he  worked  out 
his  theory  of  the  buckling  of  orthotropic  circular  cylindrical  shells. 

The  fact  that  tests  in  very  rigid  and  in  very  elastic  testing 
machines  lead  to  the  same  buckling  stress,  and  that  consequently  the 
Tsien  criterion  must  be  considered  invalid',  has  been  confirmed  for 
circular  cylindrical  shells  by  Horton,  Johnson  and  Koff  in  1961  and  by 
Almroth)  Holmes  and  Brush  in  1964.  The  same  proof  was  brough  recently 
for  complete  spherical  shells  by  Carlson,  Sendelbeck  and  Hoff  (1965). 

The  sensitivity  of  axially  compressed  circular  cylindrical  shells 
to  small  initial  deviations  from  the  exact  shape  was  demonstrated  by 
Babcock  and  Sechler  (1962,  1963)  when  they  tested'  a  series  of  very 
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accurately  fabricated  .shell  specimens  with  built-in  and  carefuly  measured 


deviations..  The  manufacturing  procedure  used  in  these  investigations  had 
originally  been  introduced  by  Thompson  in  i960  who  had  produced  thin 
spheres;  by  the  electroplating  method.  With  the  best  specimens  of  this 
v  kind  Babcock  and  Sechler  reached  J6  percent  of  the  classical  critical 
stress  when  the  a/h  ratio  was  890;  The  values  of  p~  acr/ acl 
obtained  by  Almroth,  Holmes  and  Brush  ranged  from  0.43  to  0,73.  Even 
higher  values,  up  to  0.9,  were  reported  by  Tennyson  (1963,.  1964)  when 
the  a/h  ratio  of  specimens  made  of  a  photoelastic  material  was  between 
100  and  170. 
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FIG.  7  .  Comparison  of  ,p  Values  Obtained  from  Donnell's  and  Sanders '  Equations  for  Case  SSI 
(From  International  Journal  of  Mechanical  Sciences) 


FIG.  8.  Buckles  Covering  the  Entire  Surface  of  a  Shell  Provided  with 
Close-Fitting  Mandrel  (From  International  Journal  of  Solids 
and  Structures) 
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FIG.  9*  Postbuckling  Curves  Calculated  by  Authors  Listed  in  Table  2 
(From  SUMER  No.  221). 
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FIG.  10.  Yoshimura  Buckling  Pattern 
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Load -Displacement  Diagram  for  Imperfection-Sensitive 
System  (Unstable  Branching  Point)  S  =  lOOr)3  -  200t)2+  105t) 


FIG.  15.  Load-Displacement  Diagram  for  [Imperfection-Insensitive 
System  (Stable  Branching  Point)  S  =  100^3+105J1  • 
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